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A time-accurate two-dimensional fluid code is used to compute test times in shock tubes operated at supersonic 
speeds. Unlike previous studies, this investigation resolves the finer temporal details of the shock-tube flow by 
making use of modern supercomputers and state-of-the-art computational fluid dynamic solution techniques. 
The code, besides solving the time-dependent fluid equations, also accounts for the finite rate chemistry in the 
hypersonic environment. The flowfield solutions are used to estimate relevant shock-tube parameters for laminar 
flow, such as test times, and to predict density and velocity profiles. Boundary- layer parameters such as 
3*, and 7 W , and test time parameters such as t and particle time of flight t f , are computed and compared 
with those evaluated by using Mirels’ correlations. This article then discusses in detail the effects of flow 
nonuniformities on particle time-of-flight behind the normal shock and, consequently, on the interpretation of 
shock-tube data. This article concludes that for accurate interpretation of shock-tube data, a detailed analysis 
of flowfield parameters, using a computer code such as used in this study, must be performed. 


Nomenclature 

a = speed of sound, m/s 

C p = specific heat ratio 
d, D = shock-tube diameter, m 
L - length of driven section, m 

/ = distance downstream of the shock front, 0 at 

shock front, l c at contact surface, m 
= maximum separation between the shock front and 
the contact surface, m 
m = mass flow, kg/m 3 

P = pressure, N/m -2 

T = translational temperature, K 

t f ~ particle time of flight, /'<• (fc,/p*£/,) d/. s 
V s ~ shock front velocity in the laboratory coordinates 
(wall stationary), m/s 

u = flow velocity in shock stationary frame, m/s 
W = uju, (i = 

X = nondimensional axial distance, 

u rJ) t/l c . m , t = /'. d Hu, 

y = radial coordinate, 0.0 at tube center, m 
/3 = constant 

= displacement thickness, 

Jo (up - u,p e )!(u r pX dv, m 

d t = boundary-layer thickness based on velocity, v at 
(U, - U)/(U S - «,) = 0.99. m 
d p = density displacement thickness, 

.ft (pip. - 1) m 
(a = coefficient of viscosity, kg/m s 
p - density, kg/m 3 

t = test time at / = / f , l e lu w , s 

t* = test time at / = / rm , l Cim fu w , s 
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T 

- 

nondimensional test time, r/r x 

X 

= 

hjh . 

Subscripts 


C 

= 

contact surface 

e 


downstream of the shock front 

e , 0 

= 

immediately behind the shock front 

m 

= 

limiting flow regime 

s 


shock 

St 

= 

standard conditions, 1 atm, 294 K 

u 


velocity 

w 

= 

wall 


1 = upstream of shock front 

2 = conditions between shock front and contact 

surface 

» = upstream of the shock front 

I. Introduction 

S HOCK tubes, by producing high-temperature gases, pro- 
vide a very powerful tool 1-4 for the study of gases at 
temperatures beyond those obtainable in furnaces. The shock 
tubes produce a homogeneous gas sample with enthalpy and 
pressure that can be dependably and literally calculated from 
the measured shock velocity and the conservation laws. It is 
this dependable, accurate production of gas samples that makes 
shock tubes a most attractive tool for high-temperature re- 
search. Under ideal conditions, such as a flow without wall 
or real-gas effects, the shock wave and contact surface would 
move with a constant velocity and the flow between them is 
uniform and one dimensional. However, in a real shock tube 
the homogeneity of the gas sample produced begins to suffer 
as the boundary layer grows and as real-gas effects become 
important. The interface between the driver gas and driven 
gas sample is typically turbulent. This turbulent region, along 
with boundary-layer effects and other complexities, usually 
engulfs a significant part of the heated driven gas, in some 
cases up to half of the total driven gas. When the separation 
between the contact surface and the shock front is larger than 
one tube diameter, a large portion of the heated gas sample 
is unaffected by the turbulent contact surface. 5 In these large 
tubes, this region is typically homogeneous with a planar shock 
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front. The achievement of a planar shock front followed by 
a homogeneous gas sample is essential for the validity of the 
conservation laws and the accuracy of the gas state calcula- 
tions. The interpretation of shock-tube data measured in this 
homogeneous test gas sample is straightforward, since the 
time of flight of the particles participating in the thermokinetic 
processes is accurately predictable. 

However, as the length-to-diameter (LID) ratio of the driven 
section is increased and as the initial pressure in it is reduced, 
the wall boundary-layer effects become more pronounced. 
The flow nonuniformities between the shock front and the 
contact surface add further complexities to the boundary-layer 
effects. The separation between the contact surface and the 
shock front / r increases with an increase in LiD until it reaches 
a maximum, at which point a limiting flow is established. At 
this point the shock and the contact surface move with a 
constant velocity. Duff, 5 to explain the limiting flow, referred 
to the contact surface as a leaky piston, such that a fraction 
of the test gas leaks through the wall boundary layer. He 
explained that the limiting flow occurs when the flow through 
the shock front is equal to this leakage. Roshko 7 and Hooker 8 
later studied this effect both experimentally and theoretically. 
Based on the experimental data, Roshko was able to confirm 
the validity of the basic theory. Analytical studies by Roshko, 7 
subsequently expanded and modified by Hooker, 8 are pri- 
marily based on the computations of boundary-layer param- 
eters such as , <5*, etc., in terms of relevant shock-tube 
operational parameters, such as W = u„/u ca — p<r.Jp~, the 
shock strength parameter P 2 !P\* the pressure ratio across the 
shock, and the shock Mach number. Using the boundary- 
layer parameters, the incoming mass across the shock and 
outgoing fluid mass across the contact surface can 
be estimated. In the limiting flow case when the incoming 
mass is equal to the outgoing mass across the contact surface, 
it is possible to derive a set of equations, which can be solved 
to compute the l cm and other relevant test time param- 
eters. Roshko 7 used a simplified model assuming that 1) the 
contact surface moves with a constant velocity and 2) the 
freestream has a uniform velocity profile. Later, this model 
was modified by Mirels 9-14 to account for the acceleration of 
the contact surface, flow nonuniformities, and real-gas ef- 
fects. 

Comparisons with available experimental data in all of these 
investigations were limited to shock velocities below 2.0 km/ 
s. In the early 1960s, Kemp 15 extended these models for higher 
velocities (6.0-10.0 km/s) and found that the actual measured 
test times for these velocities were, on average, one order of 
magnitude smaller than the values predicted by this theory. 
Kemp pointed out that this discrepancy may be caused by 
some additional phenomena not accounted for in these models, 
such as mixing and density gradient instabilities at the contact 
surface, and the predicted values should be regarded as the 
upper limit of the test time that can be achieved in a shock- 
tube facility. 

In these early analyses the separation between the surfaces, 
in distance or time, was taken to be as a measure of test time. 
However, the early experimentalists and developers 2 * 5 of low- 
density shock tubes did realize that the entire separation length 
is not filled with a homogeneous gas. Kantrowitz 5 thought 
that the homogeneous portion would account for only one- 
half of the separation length, or even less, as the gas tem- 
perature is increased. 

Since most of the shock-tube data relate to the kinetic pro- 
cesses in the test gas sample, the experimentalist must eval- 
uate the time required for a shock-heated particle to reach a 
specified distance downstream of the shock front for proper 
interpretation of the data. If the conditions behind the shock 
front are uniform, then the time of flight of a given particle, 
at a distance / from the shock front, is equal to 

tf = (p'MWU,) 


Here, p f „ and p are the densities immediately downstream 
and upstream of the shock front, respectively. In a nonuni- 
form flow case 



d/ 


Mirels' computations 14 show that in the limiting flow case a 
particle may take about four times longer time to reach a 
distance / = 0.9 l c m than predicted by the uniform flow case. 
Thus, in order to interpret the shock-tube data correctly, the 
experimentalist must have accurate knowledge of the flow- 
field behind the shock front. 

These findings are not new and over the years, many au- 
thors, notably Mirels, have written numerous papers to this 
subject. However, both experimental and theoretical inves- 
tigations of the homogeneity of the test gas sample are difficult 
to conduct and, for that reason, many investigators tend to 
ignore these phenomena when interpreting the shock-tube 
data. 

In this article the authors revisit the problems associated 
with inhomogeneity of the test gas in shock tubes using state- 
of-the-art computational fluid dynamic (CFD) techniques. 
Taking advantage of modern supercomputers and CFD tech- 
niques, we are able to resolve the finer details of the boundary 
layers and are able to computationally model the complexities 
of finite rate chemistry and thermal nonequilibrium. In the 
present computational study we have obtained the flowfield 
solution for a low-density operating condition with a laminar 
boundary layer, for which the assumptions made by Mirels 
are valid, namely, 1) modest real-gas effects T 2 ~~ 3500 K, 2) 
low Mach number M s ~ 8, and 3) cold wall T w - 294 K. At 
first we will compare our computational results of separation 
length (test time) with those based on Mirels’ correlations and 
experimental data of Camm et al., 4 Duff, 6 Hooker, 8 Roshko, 7 
Wuster et al., 15 and Sharma et al. 16 Then we will discuss the 
homogeneity of the test gas sample by examining in detail the 
computed flowfield profiles and its impact on the particle time 
of flight and, consequently, on the interpretation of the shock- 
tube data. Finally, we will assess the role of this new com- 
putational tool in assisting the experimentalists in interpre- 
tation of shock-tube data. 


II. Mirels 7 Correlation Formulas 
for Shock-Tube Flow 

First, we describe the basic features of boundary-layer de- 
velopments in a shock tube, although they have been pre- 
sented by Mirels and others 612 in numerous publications. The 
shock-tube flow, as discussed by Mirels, is viewed in a co- 
ordinate system in which the wall moves with a velocity 
(which equals U s in the laboratory coordinate system. Fig. 1). 
In this shock stationary system the local flow parameters be- 
tween the shock front and the contact surface are denoted by 
subscript e and the values directly behind the shock are de- 


Shock front Boundary layer 



Fig. I Shock-tube flow in shock stationary coordinates; u and U are 
the velocities in the shock stationary and laboratory coordinates, re- 
spectively. 
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noted by an additional subscript 0. The contact surface is 
located at a distance I c from the shock front. The conditions 
upstream of the shock front are denoted by subscript so and 
the conditions at the wall by subscript w. 

Duffs leaky piston analogy can be visualized by exami- 
nation of the sketch in Fig. 1. The mass flow through the 
shock front m s can be expressed by 

rh s = (p { ,u c ) () A = pJJ / 1 ( 1 ) 

where A is the cross-sectional area of the shock tube. If it is 
assumed here that the boundary layer is thin compared to the 
tube diameter, the leakage at the contact surface is given by 

m, = Lp H .„(u H . ~ u rM )S (2) 


where L is the perimeter of the tube, 8 is the characteristic 
boundary-layer thickness at l c (see Fig. 1), and p* <> and u„ - 
u cAt are the boundary-layer density and velocity, respectively. 
The added momentum due to the excess mass Am = m s — 
m ( causes the acceleration of the contact surface as observed 
in many experiments. 13 In the limiting flow (/ = l cm . Am 
0.0), the in-flow through the shock front is equal to the mass 
leakage through the contact surface. The condition rh s = m c 
provides a unique relationship between the shock-tube pa- 
rameters and the boundary-layer parameters at / — l cm and 
the maximum separation length l cm itself. 

For a laminar flow the characteristic boundary-layer thick- 
ness, 8 at / = l c m , can be written as 11 


* = P 


” U e.o)_ 


Here, the p is a constant that must be determined. Assume 
that the temperature upstream of the shock front is at a stan- 
dard condition so that T * = T st , a *. = c st , and p* = p st . It 
is further assumed that the wall remains at the initial tem- 
perature, i.e., T w = 7^. For most of the shock-tube operations 
this assumption is valid. 12 

Equations (1-3) can be combined to yield 


Pst 

p x d 2 


1 

( — 

WM S \ 

(—) 

16/3 2 

\Pc.0 

W - 1/ 

\ M / 


(4) 


however, they assume a constant freestream between the shock 
front and the contact surface. 

III. Fluid Flow Model 

As mentioned in Sec. I, the present study is based on an 
axisymmetric flowfield simulation of the shock-tube flow. A 
detailed description of this code has been presented else- 
where. 17 However, for completeness, its basic features are 
given in the following text. 

Fluid flow through the shock tube is modeled by coupling 
the thin-layer, Navier-Stokes equations with chemical reac- 
tion rate equations. The gas model includes the three major 
species (N 2 , N. and He) and accounts for finite rate chemical 
processes. A separate equation for vibrational energy is in- 
cluded, but thermal equilibrium is enforced for the present 
work. Species diffusion effects are not included. The problem 
is assumed to be axisymmetric. The contact discontinuity is 
initially assumed to be planar and the boundary layer is treated 
as laminar. The equations are solved to provide a time-ac- 
curate solution. As a result, the temporal evolution of bound- 
ary-layer growth, mass leakage through the boundary layer, 
acceleration of the contact surface, and deceleration of the 
shock front can be monitored. 

IV. Numerical Method 

The axisymmetric gasdynamic equations are solved by using 
an explicit finite volume form of the Harten-Yee upwind 
total variation diminishing (TVD) scheme, 18 which gives sec- 
ond-order spatial and temporal accuracy. The simulations 
cluster grid points at the shock and contact discontinuity and 
convect this clustered grid with these features as they travel 
down the driven tube to minimize numerical errors. 

The solutions are advanced at a Courant- Friedrichs- Lewy 
(CFL) number less than one based on the inviscid gasdy- 
namics. In the boundary layer, this CFL number can be sig- 
nificantly larger than the stability bound required by the vis- 
cous terms. To avoid the limiting time step dictated by the 
viscous terms, these terms are treated implicitly. This requires 
a block tridiagonal matrix inversion along each line of cells 
normal to the wall. The cost of this inversion is more than 
offset by the greater allowable time step. The source terms 
representing the finite rate chemical kinetics and vibrational 
relaxation are also treated implicitly. This implicit treatment 
reduces the formal time accuracy to first order. 


Here, W » uJu cM = p r . 0 /p^ and M s = UJa^ is the shock 
Mach number. 

Equation (4) can be used to compute the maximum achiev- 
able separation distance l c m , however, it requires an accurate 
estimate of p. Numerous schemes for computations of p have 
been proposed over the years 6 7 - 1113 varying in level of intri- 
cacy. Mirels, 11 using the local similarity approximation based 
on the computation of displacement thickness, has proposed 
the following expression: 


1.796 + 0.802 W \ = 

zw - 1 / ‘ ~ p„p H 


1 7,^2 _ P 2 

Z 2 T 2 * P2RT2 


(5) 


Here, C t . is the correction factor for variable pp. y Z = (y + 
\)/(y - 1) for an ideal shock. The expression for the constant 
P [Eq. (5)] includes the effects of freestream nonuniformi- 
ties. 11 Mirels recommends that for a strong nonideal shock Z 
be replaced by W (W > Z). For a shock Mach number greater 
than 8 the error caused by replacing Z with W is less than 
1%. In this study, we have used Eqs. (4) and (5) to compute 
l cm . Other correlations that were used to estimate various 
boundary-layer parameters are included in the Appendix. These 
correlations are from Ref. 12 and are valid for variable pp. 


V. Results 

Numerical results generated by this code have been com- 
pared against experimental data and were reported. 17 The 
experimental data were obtained at NASA Ames' Electric 
Arc-Driven Shock-Tube Facility. 

A. Computed Flowfield 

The shock-tube geometry used in the present study had a 
driver of 0.76 m in length and a driven tube of 12.00 m in 
length with 0.05 m i.d. The single run condition investigated 
had a driven tube fill pressure of p * = 0.525-torr nitrogen at 
294 K. The driver pressure and temperature were 6.93 atm 
and 700 K, respectively. The axial distance was divided into 
400 grid points and the radius into 63 grid points. The radial 
grid system was clustered at the wall to provide better reso- 
lution of the boundary layer. The grid in the axial direction 
was clustered near the shock front and the contact surface. 

To ascertain that the computational results are grid inde- 
pendent , a grid convergence analysis was conducted . The min- 
imum grid size at the wall was reduced in four consecutive 
runs, from 50.0 to 6.25 pm and the resultant flowfields were 
compared. All runs were started at the point in time when 
the diaphragm ruptures and carried through the point of (or 
very close to) limiting flow. The point of limiting flow was 
determined by plotting the value of l c as a function of time 
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Fig. 2 Separation distance between the shock front and the contact 
surface, / as a function of time. 


Equilibrium Riemarm Solver 
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Contact Discontinuity 1 
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Equilibrium Riemann Solver 
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Time (sec) 


i r 
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Fig. 3 Shock-front and contact surface speeds as functions of time 
in laboratory coordinates. 


and by monitoring the value l c in time. It was assumed that 
limiting flow was achieved when l c remained constant with 
further increase in time. The results of the four computations 
are shown in Fig. 2. The l cm value changed by 6.4% when 
the grid size at the wall was reduced from 50 to 25 ^tm. The 
change in l c m value was about 4.4% when the grid size was 
changed from 25 to 12 The difference between the l c m 
values for grids with 6.25 and 12.5 jum was found to be only 
2%. The displacement thickness, which is at the heart of 
Mirels' theory for boundary-layer parameters and test times 
in shock tubes, was also found to be within 2% for these two 
cases. At this point 6.25-^tm grid size at the wall was selected 
as the final grid. For the final grid, the cell Reynolds number 
at the wall as defined by the following relationship was found 
to be 3.5 at the shock front decreasing to 0.6 at the contact 
discontinuity: 

- _ pjtcviAy 

Ke C{ . i{ — 

Here, m ccU is the velocity at the center of the first cell above 
the wall and Ay is the height of the first cell above the wall. 
As a rule of thumb, Re Ci:il must be less than one to achieve 
an acceptable resolution of the boundary layer. With Re ccll = 
0.6 at the contact discontinuity, we are confident that the 
boundary layer is sufficiently resolved. However, since the 
boundary layer is thinner near the shock front, it is less re- 
solved in this region. This is revealed by the value of Re cc} , 
= 3.5 behind the shock. Consequences of this fact will be 
discussed later. To ensure that the 63 points used across the 
radius of the tube were sufficient, an additional computation 
with the 25-/Ltm wall spacing was made using 126 points across 
the radius. The results of this run are also presented in Fig. 
2. This run was not carried out as far as the other computa- 
tions, but shows little difference from the 63-point solution. 
The code requires 12.2-/xs/grid point/iteration on a Cray C- 
90. The final solution with the 6.25-ptm spacing at the wall 
was run for 760,000 iterations, requiring approximately 65 h 
of CPU time. This computation was the most expensive, since 
the time step size is controlled by the wall spacing. 

Also marked in Fig. 2 are four points corresponding to 
separation lengths of 0.5, 0.75, 0.96, and ~/ f ,„ for the finest 
grid solution. The first point is reached in 0.32 ms after dia- 
phragm rupture, the second in 0.66 ms, the third in 2.01 ms, 
and the fourth at 3.35 ms. These four points will be referred 
to later. Mirels’ correlations predict / t m = 0.152 m, about 
0.8% higher than the value of 0.1508 m based on our final 
solution (t = 3.35 ms). It is conceivable that if the code were 
allowed to run further, the computed value of l cm may reach 
Mirels’ value, however, it would more than double the in- 
tegration period and was deemed impractical for this study. 

The shock front and contact surface velocities (in the lab- 
oratory coordinates) as a function of the time from diaphragm 
rupture are shown in Fig. 3. The deceleration of the shock 



Fig. 4 Translational temperature as functions of axial distance in 
laboratory coordinates. 


front and acceleration of the contact surface are clearly seen 
in this figure. Because of the large deceleration of the shock 
front, the particles passing through the shock in the early 
stages of the flow are processed by a shock whose strength is 
significantly greater than its strength at later times. This phe- 
nomenon causes additional inhomogeneity in the test gas sam- 
ple, besides those due to the boundary-layer effects discussed 
earlier. 

The temperature profiles over the whole shock tube are 
shown in Fig. 4. Plots for the four points indicated in Fig. 2 
are shown in Fig. 4. For clarity, the plots on Fig. 4 are shown 
in the laboratory coordinates. Two observations are made 
here: 

1) The difference between the temperatures directly down- 

stream of the shock front and at the contact surface is much 
larger at point no. 1 (/,. = 0.5/ ( A T ** 458 K) than that at 

point no. 4 (/,. = A T = 153 K). 

2) A large temperature gradient exists upstream of the con- 
tact surface at point no. 1, which as the shock moves further 
along the driven tube, almost disappears at l c = The 
translational temperature at the contact surface at point no. 
1 is about 4300 K, as compared to 3842 K at the shock front. 
Fluid particles at the contact surface were processed before 
those at the shock front, and as seen from Fig. 3, therefore 
processed by a stronger shock. The temperature of 4300 K at 
the contact surface at point no. 1 is consistent with the post- 
shock temperature for the initial shock speed of 3200 m/s. As 
the gas sample moves along the shock tube some of the fluid 
particles at the contact surface leak out through the boundary 
layer, resulting in a slightly cooler contact surface at later 
times (see point no. 4 in Fig. 4). 

Velocity contours corresponding to the limiting flow solu- 
tion are shown in Fig. 5, which clearly shows the shock front, 
contact surface, and the boundary-layer buildup. Each con- 
tour in this figure represents a 30-m/s change in velocity. The 
flow streamlines in the shock fixed coordinates are shown in 
Fig. 6. It is very instructive to see that, this being the limiting 
flow, all of the streamlines (and the mass flow through them) 
at the contact surface flow into the confines of the boundary 
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Fig. 5 Velocity contours in the test gas sample; P a = 0.525 torr, T„ 
= 294 K, and d — 0.05 m. Driver: T — 700 K and P — 6.93 atm. 
Shock stationary coordinates. 



Fig. 6 Streamline contours in the test gas sample; P 0 = 0.525 torr, 
T m = 294 K, and d = 0.05 m. Driver: T = 700 K and P = 6.93 atm. 
Shock stationary coordinates. 
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Fig. 7 Boundary-layer thickness parameter as a function of distance 
measured from the shock front. 


layer. This clearly demonstrates the leaky piston analogy of 
Duff, as was mentioned earlier. Since no experimental veri- 
fication of the leaky piston concept has been produced, it is 
instructive to visualize and demonstrate leaky piston phenom- 
ena using the present code, supporting the theories of Mirels, 
Duff, and Roshko. 

B. Mirels’ Theory vs Present Analysis 

One of the most interesting relationships that appeared in 
Mirels’ work was the universal relationship between the non- 
dimensional maximum shock/contact surface separation dis- 
tance and the shock Mach number. Mirels' universal rela- 
tionship is based on several other parameters, such as, S u , 
8*, 8 P , t k , q H ., and other shock-tube parameters, such as t 
and X along the driven length. The expressions for the bound- 
ary-layer parameters are listed in the Appendix. The test time 
parameters are defined as follows: 



The S u for the limiting flow, based on Mirels’ correlation as 
well as from the present study, is plotted in Fig. 7. The value 
based on Mirels’ correlation is constant over the distance since 
the correlations assume a uniform freestream ( u e = w^o)- 
However, since in the present solution the freestream con- 
ditions are changing as a function of distance, the computed 
boundary-layer parameters do not remain constant over the 
distance l c between the shock front and the contact surface. 
Ideally, the values of these parameters from both sets should 
be the same at l — 0.0, immediately behind the shock front. 
However, at / = 0.0 the boundary layer is very thin and, even 
with the grid size of 6.25 /im at the wall, it is not possible to 
resolve the boundary-layer parameter near the shock front 
accurately. Also, the shock front tends to be curved near the 
wall and the resultant freestream near the wall has its own 
effect on the boundary-layer growth. For these reasons, the 
values of the boundary-layer parameters do not agree with 


Table 1 Values for / = 0.9 l cjm at limiting 
flow regime" 



Present 

study 

Mirels' 

data 

Difference, 

% 

S„ 

4.04e-4 

2.85e-4 

41.7 

8* 

8.74e-4 

7.03e - 4 

24.3 

K 

9.03e - 5 

7.87e -5 

14.7 


2.02e + 3 

2.00e + 3 

1.0 


6.02e + 6 

6.27e + 6 

4.0 


"Here. S„, 5*, and arc in the units of m i : . f is in the 
units of N m"*-, and q n is in the units of W m - v: . 



Fig. 8 Nondimensional test time as computed and based on Mirels' 
correlations compared with experimental data; U t < 3.0 km/s. 



Fig. 9 Nondimensional test time as computed and based on Mirels’ 
correlations compared with experimental data; U s > 3.0 km/s. 


Roshko’s and Hooker’s data show reasonable agreement with 
the theoretical curves, especially the asymptotic nature of f. 
Roshko used the value of /3 = V3 in analyzing his data. The 
the values based on Mirels’ correlations at / = 0.0. If the 
results for / > 0.08 m, where the boundary layer is well re- 
solved, are extrapolated (see Fig. 7) to the shock front, then 
the agreement is good. It should be noted that the effect of 
nonuniform freestream on 8 U is particularly significant (see 
Table 1). 

For computing t, l cm was computed using Eqs. (4) and 
(5). The computed value of l c m is 0.1508 m against 0.1519 m 
based on Mirels’ correlations (a difference of 0.7%). As men- 
tioned earlier, Mirels’ correlations for l cm are corrected for 
the nonuniformities in the freestream and the correlations for 
8 * , etc. , are not. That is why, in spite of a significant difference 
in the values of the boundary-layer parameters, especially 5*, 
there is an excellent agreement in the values of l c It is not 
surprising that a plot of f vs X using Mirels’ theory also shows 
an excellent agreement with the present computations (Fig. 
8). Also shown on the plots are the experimental data of 
Duff, 6 Roshko, 7 and Hooker 11 reproduced from Ref. 11. 
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parameter f is very sensitive to the value of /3 near its asymp- 
totic value (f — 1.0). If & were correctly evaluated Roshko’s 
data should be asymptotic, i.e., t — > 1.0 (since l cm =< /3 -2 ). 

It should be noted that all of the experimental data plotted 
in Fig. 8 are for M < 9 and the present computations are for 
M = 8.2. 

As mentioned earlier, at higher velocities the flowfields are 
affected by some additional phenomena 13 not accounted for 
in the present code or in Mirels’ correlations, such as mixing 
and thermal instabilities. As a result the recorded test times 
are much smaller than those predicted by the theory. Some- 
times these instabilities cause the disintegration of the contact 
surface and under these conditions it is difficult to make any 
reliable measurement of the test time. A plot of test times at 
higher velocities from the experimental data of Sharma et 
al., 16 Wuster et al., 15 and Camm et al. 4 is shown in Fig. 9 
along with the theoretical predictions. For computing r , = 
/ cm /w H „ l cm for all of the experimental data plotted in Fig. 9 
was taken from Mirels’ computations (Ref. 11, Fig. 6). In the 
present simulation, for simplicity, assume a clean “breakup” 
of the diaphragm, which is rarely the case. Also, it is difficult 
to simulate the shot-to-shot variations caused by minute im- 
perfections in driver assembly, such as diaphragm defects, 
trigger wire geometry, fill gas pressures, etc. Experience 16 
shows that these unavoidable imperfections affect the way the 
diaphragm opens and in certain cases cause the onset of var- 
ious interface instabilities. These phenomena are reflected in 
the test data presented in Figs. 8 and 9. In view of these facts 
one concludes that the theoretical curves provide only the 
upper limit for the test time, and in real shock tubes the test 
time is usually smaller. 

In summarizing the results for laminar flow with relatively^ 
small real-gas affects, Mirels’ universal correlations provide 
good estimates of the test times. The accuracy of these pre- 
dictions is comparable with that achieved by using an intricate 
tool such as the present fluid code. However, the boundary- 
layer thickness is very sensitive to the flow nonuniformities 
and only a computer code such as is presented here or detailed 
experimental measurements can provide the needed infor- 
mation about it. 

C. Test Gas Homogeneity 

Figure 10 depicts the relative variation of various ther- 
modynamic parameters, namely, 7\ P, and u e , with respect 
to their magnitude at the shock front (7^, P s , and u e 0 ), as a 
function of normalized distance from the shock front, lll cm . 
The properties shown correspond to the limiting flow regime 
(/ = l cm ). The flow velocity in the shock stationary coordi- 
nates reduces to 0.0 as a particle moves towards the contact 
surface, which is as expected in the limiting flow case. The 
changes in temperature and pressure are less dramatic. Pres- 
sure increases by 9.5% and the temperature by 4.8%. Figure 
11 shows the variation of these properties at a location when 
the separation between the contact surface and shock front 



at the limiting flow. 



I /(0.5* W 


Fig. 11 Normalized 3\ P, and u, as functions of normalized distance 
at a location where / = (0.5/,. m ). 


is about l cm l 2 ( t = 0.32 ms). Even at this early stage, the 
reduction in u e can be as high as 60% . 

As was discussed in Sec. V.A, the temperature profile in 
the test gas sample is a very strong function of the shock speed 
history. At the early stages of the flow the particles are sub- 
jected to a large temperature gradient (Figs. 2 and 3), and 
consequently, the temperature variations in the gas sample 
in Fig. 11 are mucl^larger than in the limiting flow shown in 
Fig. 10. 

D. Particle Time of Flight 

The inhomogeneity of the test gas sample (seen in Figs. 10 
and 11) affects the particle time of flight. As a result, for 
proper interpretation of shock-tube data, one requires the 
pertinent information about the flowfield in the test gas sam- 
ple so that the necessary corrections to the data can be made. 
For a uniform flow downstream of the shock front, the particle 
time of flight is given by t f = (//l/,)(p, <0 /p*) = t f 0 . In the case 
of an inhomogeneous gas sample, it is equal to 



For strong shocks, Mirels 14 suggests the following correlation 
for the computations of 


f = l = 

tfO Pe.O 1 



Mi - (///c.J 1/2 ]- ! 




(7) 


To visualize the effect of inhomogeneity on the particle time 
of flight, the nondimensional parameter t f !t f0 is plotted in Fig. 
12 as a function of distance from the shock front. This plot 
is produced for the limiting flow regime using the data for 
the centerline of the tube. Computations using Eq. (7) are 
also plotted in the figure. For the plot based on the present 
computations the velocity data u f correspond to the centerline 
grid point. There is a fair agreement between the present 
computations and the plot based on Mirels’ correlation, for 
Ul cjn < 0.6. For lll c m > 0.6, the data from Mirels’ correlation 
for the particle time of flight are on average about 15-20% 
higher than the present real-gas solution. At Hl cm = 0.9, the 
values for t f /t f0 based on the present computations and on 
Mirels’ correlations are 4.2 and 4.5, respectively. In other 
words, a particle located at the centerline will take about 4.2- 
4.5 times longer, due to varying thermodynamic conditions 
in the gas sample, than in a uniform gas sample (if u e were 
constant from the shock front to the contact surface). 

As was mentioned, the plots in Fig. 12 are for the limiting 
flow regime. Computations for a location corresponding to 
point no. 1 on Fig. 4 (/ = 0.5/ c/ „) indicate that even at this 
early stage the ratio t f !t ^ has a value larger than 1.00 and 
gradually increases with increase in distance as measured from 
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Fig. 12 tfltfo as a function of normalized distance at limiting flow. 



Fig. 13 as a function of distance along the shock tube. 

the shock front. At this location the particle takes about 1.8 
times longer than in a uniform flow case to travel a distance 
equal to 0.5l cm . 

In Fig. 13, the evolution of the nondimensional parameter 
(tfltfo) c with time can be seen. Here the subscript c denotes 
the fact that the times of flight are for the distance l c . The 
parameter ( t f /t f0 ) c increases almost linearly as the shock front 
moves along the shock tube. The maximum value of 
is about 100. It must be remembered that strictly speaking, 
the flow solution has not achieved the true limiting flow con- 
dition. It is very near to it and will still require several mil- 
liseconds to achieve it. Ultimately, when the limiting flow 
regime is achieved, the parameter (t f /t f0 ) c should become in- 
finite, as Mirels’ correlation predicts. 


1) The experimentalist needs to know the thickness and 
thermal properties of the boundary layer through which the 
laser beam passes or through which emission measurements 
are being made. 

2) In searching through operating conditions of the shock 
tube or for an expansion tube, flowfield information behind 
the contact surface may be required. 

3) Flowfield properties of the inhomogeneous gas sample 
may be required for proper evaluation of particle time of 
flight. For all of this information the flowfield solution of the 
shock tube in question will be essential. 

Mirels’ universal correlations are excellent tools for esti- 
mating the test times in shock tubes, especially for the limiting 
flow regime. However, as we have shown, for correct as- 
sessment of the boundary-layer parameters, one requires a 
detailed flowfield solution of the shock tube. With our present 
knowledge of real-gas effects, modern CFD techniques, and 
the availability of supercomputers, it is now possible to pro- 
vide improved data to experimentalists for proper interpre- 
tation of shock-tube data. 

VII, Conclusions 

Capabilities of a time-accurate shock-tube flow code have 
been demonstrated. The boundary-layer parameters, such as 
5 U , 5*, S p , f K (see Appendix), and other shock-tube pa- 
rameters such as t, have been computed using this code. The 
parameters such as rand l cm agree well with the values com- 
puted using Mirels’ correlations. However, the boundary-layer 
thicknesses 5 1( , and 5*, do not. The particle time of flight 
data for the limiting case also show good agreement. 

The computations also indicate that flow nonuniformities 
do exist even at early stages of the shock-tube flow, and 
consequently, one must make correction for particle time of 
flight in data taken at those locations. However, to do so one 
needs to make use of computer codes such as the one pre- 
sented in this article, since Mirels’ work does not provide 
correlations for cases other than the limiting flow case. Most 
shock tubes in use are not long enough so that the limiting 
flow conditions can be achieved. Therefore, in most cases, it 
is not possible to use Mirels' correlations. 

Mirels’ correlations and the present shock-tube code have 
been validated against each other. At the same time, it is 
concluded that for proper interpretation of shock-tube data, 
one must make use of a code like the one presented in this 
article. 


VI. Discussion 

Many important sets of experimental data 19-33 have been 
acquired in shock tubes and these facilities continue to play 
a very important role in investigations of chemical kinetics, 
aerothermodynamics, and space technology. The experimen- 
talists working with shock tubes have been continually search- 
ing for improved and accurate techniques for data interpre- 
tation. Therefore, whenever new data interpretation tools are 
available, it is imperative that we assess their capabilities. 

In line with this philosophy, the main objective of this study 
is to demonstrate the usefulness of the present computer code 
in analyzing the flowfield in shock tubes, providing necessary 
information about the nonuniformities of the test gas sample 
to assist the experimentalist in interpretation of shock-tube 
data. The input conditions for this case are selected so that 
the basic assumptions made by Mirels also apply to the present 
flowfield. Using this solution we proceeded to validate this 
computer code against Mirels' correlations and vice versa. To 
the best of the authors’ knowledge, no other multidimensional 
CFD solution has been compared with the results obtained 
from Mirels’ correlations. 

In most shock-tube experiments the task of data interpre- 
tation requires more information about the flowfield than 
provided by Mirels’ correlations. 


Appendix: Mirels’ Correlational Formulas 

For completeness, Mirels’ correlational formulas 12 for lam- 
inar boundary layer that have been used in this study are 
reproduced. 

Boundary-layer parameters: 



The boundary-layer parameters with real-gas effects and 
corrections for variable fxp as published by Mirels 12 are given: 
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where 


c - PrfXr = — L.th z t k 

PhM^ z 2 7", Pi’ 22 P 2 /?7' 2 

is the correction factor for variable and \ = The 

parameters <£"(0), tj,,, /5 (1 - G) dr?, etc., are computed by 
the following relationships: 


0.7979 

4>\ 0) = — — (1 + 1.285c + 0.3827c 2 ) 1/4 (0.1%) 

(A6) 
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0.7979s 1 * 
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2.575c 1 * 

(1 + 0.803c + 0.1069c 2 ) 1 ' 
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where e = W - 1 and cr is the Prandtl number. The maximum 
error is indicated for each expression. 
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